Damping by branching: a bioinspiration from trees 
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Man-made slender structures are known to be sensitive to high levels of vibration, due to their 
flexibility, which often cause irreversible damage. In nature, trees repeatedly endure large amplitudes 
of motion, mostly caused by strong climatic events, yet with minor or no damage in most cases. A 
new damping mechanism inspired by the architecture of trees is here identified and characterized in 
the simplest tree-like structure, a Y-shape branched structure. Through analytical and numerical 
analyses of a simple two-degree-of-freedom model, branching is shown to be the key ingredient 
in this protective mechanism that we call damping-by-branching. It originates in the geometrical 
nonlinearities so that it is specifically efficient to damp out large amplitudes of motion. A more 
realistic model, using ffexible beam approximation, shows that the mechanism is robust. Finally, two 
bioinspired architectures are analyzed, showing significant levels of damping achieved via branching 
with typically 30% of the energy being dissipated in one oscillation. This concept of damping-by- 
branching is of simple practical use in the design of slender flexible structures. 
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I. INTRODUCTION 

Vibrations in man-made structures are a central prob- 
lem in mechanical engineering (Den Hartog 2007). They 
may result from external excitations such as wind, im- 
pacts, or earthquakes; or from internal excitations, such 
as a flow or moving parts. Their consequences are numer- 
ous in terms of functionality losses due to wear, fatigue, 
or noise, to cite a few. We distinguish here between low 
and high levels of vibration. The former, where displace- 
ments are small in comparison with the characteristic size 
of the structure, may induce some of the long-term above- 
cited consequences. The latter generally cause short- 
term failures and irreversible damage to the structure 
by fracture or plastic deformation (Collins 1993). These 
large amplitudes of vibrations may be particularly ex- 
pected in slender structures, or assemblages of them, due 
to their high flexibility. 

In the most general framework of vibration analysis, 
the amplitude of motion results, on one hand, from the 
characteristics of the loading, and on the other hand, 
from the characteristics of the structure in terms of in- 
ertia, stiffness, and damping (Humar 2002). Damping 
here refers to the capability of the structure to dissi- 
pate mechanical energy, whatever the physical mech- 
anism involved (viscoelasticity, friction or interaction 
with a fluid). A high level of damping in a struc- 
ture is a standard way to avoid large-amplitude mo- 
tions. This is generally achieved with passive techniques, 
such as the classical addition of tuned mass-damper sys- 
tems (Den Hartog 2007, p. 119), or with active or 
semi-active means such as piezoelectric materials, mag- 
netorheological fluids, shape memory alloys, or even sim- 
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pie hydraulic actuators in feedback or feedforward sys- 
tems (Preumont 2002). All these approaches have limits 
in terms of their range of acceptable deformations, dis- 
placements, or simply in terms of cost or maintenance. 
But, more fundamentally, they stem from an ad-hoc cor- 
rective point of view, rather than from a consistent design 
perspective. New approaches are needed, particularly for 
very slender and light structures, such as antennas, which 
may encounter large flow-induced amplitudes of vibra- 
tion (Pai'doussis, Price and de Langre 2011). 

Nature may give insights into highly ef&cient mechani- 
cal solutions in vibration problems, for instance in shock- 
absorbing devices (Yoon and Park 2011). Interestingly, 
slender structures are ubiquitous in nature, particularly 
in plants. Most of these plants are regularly submit- 
ted to natural flow excitations by wind or current caus- 
ing vibrations (de Langre 2008). Although these vi- 
brations contribute to some biological functions such 
as in seed or pollen dispersion, extreme events such as 
storms may cause dangerous large amplitudes of mo- 
tion (Niklas 1992). Therefore, in areas where intense 
flows are common, plants are likely to possess efficient 
damping strategies. 

From a biomimetic point of view, the dynamical be- 
haviour of trees, which has been extensively studied, 
is certainly a possible source of inspiration. Scannell 
(1984) hinted that trees might possess a "qualitative 
mechanical design principle [...] beneficial to the tree's 
survivability in conditions of strong atmospheric tur- 
bulence". Niklas (1992, p. 183) noted that "exper- 
iments indicate that branching [...] dampens natural 
frequencies of vibration". At this point it is necessary 
to clarify what is generally agreed to cause damping 
in trees. Firstly, the constitutive material, wood, is 
known to have a inherent viscoelastic behaviour caus- 
ing dissipation, this itself has been the source of bioin- 
spired material (Spatz, Emanns and Speck 2004). Sec- 
ondly the aeroelastic interaction with the surrounding 
air causes forces in the opposite direction to the local ve- 
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locity in the tree, thereby causing a strongly amphtude- 
dependent dissipation (Blevins 1990). FinaUy, when con- 
sidering the overah motion of the tree by bending of the 
trunk, another mechanism is often described as "struc- 
tural damping" (Briichert, Speck and Spatz 2003, Speck 
and Spatz 2004, James, Haritos and Ades 2006, Moore 
and Maguire 2008). This third mechanism refers to the 
possible transfer of mechanical energy from the trunk to 
the branches, where it will be eventually dissipated by the 
two aforementioned aeroelastic and viscoelastic damp- 
ing mechanisms (Sellier and Fourcaud 2009). Spatz, 
Bruchert and Pfisterer (2007) have suggested that the 
frequencies tuning of the branches plays a key role in 
this hypothesised energy transfer mechanism in trees. In 
fact, Rodriguez et al. ( 2008) analyzed the architecture 
of a walnut tree, figure 1(a) and have shown that the 
modal frequencies are close and that the modal shapes 
are strongly localized in the architecture. The former 
characteristic is classically favorable to modal energy ex- 
changes in dynamical structures and the latter is con- 
sistent with localized energy transfers from the trunk to 
the branches. In order to develop strategies for bioin- 
spired designs of slender structures including an efficient 
damping effect for large amplitudes, it is crucial to clarify 
the mechanism involved in the energy transfer that many 
authors invoke. Note that we do not intend to explore 
here the idea of the classical tuned-mass damper, which 
is based on a purely linear mechanism which, by defini- 
tion, is not amplitude-dependent. The aim of the present 
paper is therefore to identify and characterize the elemen- 
tary mechanism causing nonlinear modal energy transfer 
and damping in a branched structure. 

For this purpose, we first consider the simplest model 
of a branched dynamical system in Section |TT1 a spring- 
mass model of a Y-shape. Section IIIII shows, using a 
beam-finite-element model, that the main results of the 
previous section are also valid for a more realistic con- 
tinuous structure of a Y-shape. Based on these results. 



two illustrative designs of bioinspired slender structures 
exhibiting efficient damping-by-branching are proposed 
in Section ITVl The generality and possible extensions of 
our approach are discussed in Section [V] 



II. LUMPED-PARAMETER MODEL OF A 
Y-SHAPE 

In order to reduce the dynamics of a branched struc- 
ture to its simplest possible features, we treat the case of 
a spring-mass model of a Y-shape consisting of a trunk 
and two branches. Since we are interested in the branch- 
ing effect, viscous damping is introduced in the branches 
only. The equations of motion are written with dimen- 
sionless variables and the dynamics is studied with an 
emphasis on the damping of the whole structure. 



A. Model 

The model consists of three massless rigid bars linked 
by rotational springs and supporting three masses, fig- 
1(b) I The first bar, mimicking a trunk of length Zi, 



ure 

is linked to the ground by a rotational spring fci and 
supports a mass mi. The branches are two symmetrical 
bars of length I2, each forming an angle (ph with respect 
to the trunk axis. Each branch is linked to the tip of the 
trunk by a rotational spring ^2, and supports a mass m2. 
The motion of the trunk is defined by the angle 9, and 
we consider only the symmetrical motion of the branches 
defined by the angle 0, figure 



tion is made in order to simp 



1 (c) Note that this restric- 



ify the following dynamical 



analysis but has a small impact on the damping mech- 
anism described in this section, as it will be seen in the 
following sections for much more complex models which 
have no such restrictions. The kinetic energy is the sum 
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of the kinetic energy of each mass, 
1 



^=2 



{mih- + 2m2(/i + 2I1I2 cos((/«b + </>)+ h )) ' 



2TO2^ 



2 12 



(1) 



The potential energy is the sum of the potential energy 
of each spring, 



V^l {kit 



(2) 



The equations of motion are derived using T and V in 
the classical framework of Lagrangian dynamics (Humar 
2002). They read 



JeO + kiO = 4m2/i^2 



) sm(0b 



2m2^2^0 + 2fc2(?!) = ~2m2hl2()^ sin(0b 
where 



(3) 



Je = mi/i^ + 2m2(/i^ + 2hl2COS(j3b + h^) (4) 

and J'g{(t>) — cos(0b + (f>) ^ cos(/)b- The left-hand side 
of this system of equations represents two simple linear 
harmonic oscillators. Denoting the generalised displace- 
ment vector [0,0], the two corresponding normal modes 
of the system are directly [1, 0] and [0, 1] since there is no 
linear coupling between and cj). The two modal angular 
frequencies are respectively 



fci 
Je' 



and Lij2^ 



2k2 



2777,2^2 



(5) 



The first mode consists of motion involving 6 only, and 
the second mode involving (j) only. Therefore, in the fol- 
lowing, they are referred to as the trunk mode and the 
branch mode respectively. These two modes are coupled 
by the nonlinear terms of the right-hand side of [3J repre- 
senting the geometric nonlinearities. 

A dimensional analysis reveals the existence of four 
dimensionless parameters describing the dynamics of the 
model. We chose the dimensionless time r cjii, the 
branching angle , the ratio of angular frequencies f2 = 
ijj2liO\i and r, the ratio between the inertial terms of 
the branch mode and the trunk mode, multiplied by the 
length ratio Z1/Z2, 



r = 



2m2^2 l\ 2m2l\l2 



Je 



(6) 



The dynamics is described by the variables 6(t) = 
(){t)\/h/h and $(t) = (f>{t)- As mentioned earlier, we in- 
troduce energy dissipation in the form of a viscous damp- 
ing rate ^b in the branch mode only. The dimensionless 
equations of motion are 



e + e = 2r e$ sin((/.b + *) - e 

$ + 2f76,$ + = -6)2 sin((/)b + 



(7) 



Note that J^{(j)) = since $(r) = (f>{t). The dimen- 

sionless total mechanical energy is 



E{r) 



(2r (cos( 



COS 0b ) 



1)92 



(8) 



Since the two modes are coupled by nonlinear terms, en- 
ergy can be exchanged between them. In this case, the 
dissipation in the branch mode may damp the energy 
received from the trunk mode, resulting in an effective 
damping of the whole structure. 



B. Damping criterion 

In the following, we examine the free vibrations follow- 
ing an initial condition 



6(0), 6(0), $(0), 6(0) =[60,0,0,0 



(9) 



such that the energy is located in the undamped trunk 
mode only. This will allow us to easily demonstrate 
damping by nonlinear modal energy transfer, if any, in 
a purely linear framework, energy would remain in the 
undamped trunk mode with no way of being dissipated. 
The amplitude of the initial condition, 60, determines 
the initial energy E{Q) = Eq, using [H For the sake of 
clarity, the energy E is normalized so that the initial en- 
ergy Eq is 1 when 60 = 7r/2 corresponding to a horizontal 
trunk initial condition. Note that ground interaction is 
here neglected. 

During free oscillations, a part of the energy trans- 
ferred from the trunk mode to the branch mode is dis- 
sipated. The total energy decay over the first period of 
the trunk mode is AE = Eq — E{2'k) so that the effective 
damping rate of the whole structure can be defined as 



eff 



1 AE 







(10) 



The effective damping rate, ^off, is commonly related to 
the quality factor by Qeff = l/(2^cff)- Note that ^eff 
represents the dissipation of the whole structure and not 
that of the trunk mode. In fact, studying exclusively 
the trunk mode damping is not appropriate since energy 
transfer can be reciprocal from the branch mode to the 
trunk mode as well, as will be seen in figure [S] The total 
energy decay AE is given by the work of the damping 
term of the branch mode equation over one period of the 
trunk mode: 



AE 



2rj^b$ dr. 



(11) 



Here, the coefficient S/tt^ comes from the normalization 
chosen for E. We analyze now the effect of the initial 
energy Eq and the design parameters 0b j Cbj and F 
on the effective damping, ^off- 
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C. Energy transfer by internal resonance 

In this section, we consider a low initial energy level 
so that Qq = e, where e ^ 1 is a small parame- 
ter. The harmonic balance method (Nayfeh, Mook and 
Corporation 1979) is used with the angles and $ de- 
veloped as power series of e, 

e(r) = eQiir) + e^e2{r) + ■ ■ ■ , (12) 
$(r) = e$i(r)+e2$2(r) + ... . (13) 

The initial condition [9] requires that 

ei(0) = l, and 62(0) = $i(0) = $2(0) = 0. (14) 

Substituting [T2] and [T3l in the dynamical equations [71 and 
using [TU the first-order terms are 

9i=cosT, and $1 = 0. (15) 

The second-order terms satisfy respectively 

82-0, and $2-f21]^b$2 + f^^$2 = -ei%in(/)b. (16) 
Therefore, for small angles, [7] reduces to 

n — cos 2t) 

<^ + 2n^x,<^ + n'^<i> = -Qlsm(f)],^ — ^. (17) 

This is the equation of a simple harmonic damped oscil- 
lator, driven by an harmonic force, that can be analyti- 
cally solved (Humar 2002). A resonance exists at $7 = 2; 
since fl is the frequency ratio of the two modes, this is 
classically referred to as a 1:2 internal resonance (Nayfeh 
et al. 1979). In the following, we will discuss the infiu- 
ence of Q near this particular value. A general result for 
a forced damped oscillator is that the amplitude of mo- 
tion is proportional to the amplitude of the driving force. 
As can be seen in[T71 the amplitude of the driving force 
is proportional to 69 sin 0b, and therefore to £^0 sin 0b- 
The efi'ective damping ^cff , defined by [TU] and [TTl can 
therefore be simply expressed as 

fcff = £^orsin2 0b e(eb,0). (18) 

Remarkably, shows that ^eff increases linearly with 
the initial energy Eq: such a nonlinear damping propor- 
tional to the energy is typical of an oscillator following the 
generic equation: 6 -I- kQ^ -1-8 = (Nayfeh et al. 1979). 
Besides, [TS] shows that ^cff is proportional to sin^ 0b so 
that the effective damping is maximal for a branching 
angle 0b = 7r/2, corresponding to a T-shaped structure. 
Conversely, for a non-branched structure, where 0b = 
or TT, the effective damping is zero. The effective damping 
is also proportional to the relative modal mass ratio F. 
The dependency of f off on and ^b is shown in figure [2] 
as a contour-map of the normalised effective damping 
^, computed using a mathematic symbolic software to 
solve [TTl for $, and then using successively [TTl [TUl and 
[T5I for ^. We observe that a significant level of damping 




FIG. 2: Normalized effective damping 1181 brought by 
branching, as a function of the branch mode damping ^b, and 
the branch/trunk modal frequency ratio il. A high damping 
is found near the 1:2 internal resonance, i.e. at 17 = 2. 

is present over a wide range of parameter values. The 
effective damping shows a maximum for branch damping 
near 0.2 and frequency ratio near 2. Accordingly, the 
values 0b = tt/2, ^b = 0.2, f7 = 2 and F = 0.2 wiU be 
used as a reference in the remainder of this paper. 

As expected, we observe in figure [2] that there is no ef- 
fective damping for ^b = since mechanical energy can- 
not be dissipated in the structure. Interestingly, for any 
arbitrary small value of ^b, the effective damping is fi- 
nite. In a purely linear framework, the effective damping 
would be zero for any value of ^b, since the total en- 
ergy would be confined to the trunk mode, without any 
possible transfer to the branch mode where dissipation 
occurs. In other words, the effective damping is due to 
the geometric nonlinearities. 

D. Effects of the design parameters 

In order to obtain the full dynamics and the corre- 
sponding effective damping at any energy level with an 
emphasis on the effects of the design parameters 0b, Cb, 

and F, the dynamical system [7| is now solved numeri- 
cally by means of an fourth-order explicit Runge-Kutta 
temporal scheme. As a typical example, figure [3| shows, 
for an initial energy level Eq = 1, the evolution of the 
total energy E and of the modal energies 

= + (19) 

Eq, = -^r (^$2 {^2^2^ . (20) 

Note that the total energy E, [8l is the sum of Eq, E^ 
and a nonlinear energy term. The energy exchange be- 
tween the two modes is clearly shown. Since energy is 
dissipated in the branch mode, the total energy decays 
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FIG. 3: Typical evolution of the total energy, E ( — ), and 

modal energies, Eq ( ) and E^ (• • • ) with the respective 

evolution of the trunk angle, ( ), and branch angle, $ 

(■•■), of the spring-mass model of a Y-shape, as a function 
of time over three periods of the trunk mode. The initial 
energy is Eq = 1, i.e Go — it/2, in the trunk mode only. 
The total energy decreases as a consequence of the energy 
nonlinearly transferred to the damped branch mode. The 
design parameters are set to <^b ~ 'r/2, ~ 0.2, (1 = 2, and 
r = 0.2. 



at an effective damping rate ^eff- Figure S] shows the Eq- 
dependency of this effective damping ^off, in comparison 
with the analytical prediction of the previous section. 
As expected, the analytical approach corresponds to the 
limit of the numerical solution as Eq tends to zero. As 
Eq increases, the analytical approach increasingly over- 
estimates the numerical effective damping. However, the 
ratio ^bs/Eq — constant in the analytical approach — 
remains finite: this constitutes the essential effect of 
branching on damping. 

The influences of the design parameters 4>h, £,h, ^, and 
r on the effective damping, scaled by the initial energy, 
^es/Eo, are represented in figures [5l[a)-(d) for three ini- 
tial energy levels. As expected, the analytical and nu- 
merical approaches yield identical results for low energy 
levels and therefore are represented by the same curve 
for Eq ~ 0.01. Consistently with figure Hj we observe 
that the analytical approach overestimates the numer- 
ical effective damping as Eq increases. In figure 5(a) 



at low energy levels, the effective damping is propor- 
tional to sin^ 0b as predicted by the analytical approach 
[T51 The optimal branching angle 0b shifts from 7r/2 to 
slightly higher values when the initial energy level in- 
creases. Therefore, in order to obtain an optimum effec- 
tive damping at any energy, a good compromise would 
be to set the branching angle 4>h between 7r/2 and 27r/3. 



Ccff 2 




FIG. 4: Effect of the initial energy level, Eq, on the effec- 
tive damping, ^efi, of the spring-mass model of a Y-shape: 
( ) analytical effective damping from the low energy ap- 
proximation, 1181 ( — ) numerical effective damping from the 
full dynamics integration. The design parameters are set to 
the same values as in figure [S] 




FIG. 5: Effects of the design parameters on the effective 
damping scaled by the initial energy, 5eff/i?o, of the spring- 
mass model of a Y-shape: Eq — 0.01 ( — ); -Eo =0.1 ( ); 

and Eo = 1 (■■■)• Unless varied, the design parameters values 
are: cj)b = 7r/2, = 0.2, = 2, and T = 0.2. (a) Effect of the 
branching angle 0b. (b) Effect of the branch mode damping 
^b. (c) Effect of the branch/trunk modal frequency ratio SI. 
(d) Effect of the branch/trunk modal mass ratio F. 



Figure 



5(b) shows that a significant effective damping 



is created by a large range of branch mode damping, ^b, 
as was also found in the low-energy analytical approach. 
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Analogously, the typical shape of the Q-dependency at 
low energy is also conserved when the energy increases, 
figure 5(c) The simple F-dependency on the effective 
damping is shown in figure |5(d)[ the modal mass ratio 
r has to be maximal in order to get the highest possible 
effective damping at any energy level. 



III. FINITE-ELEMENT MODEL OF A Y-SHAPE 

The damping-by-branching mechanism described in 
the preceding section is now analysed in the case of a 
more realistic continuous beam structure of a Y-shape. 
The same approach is used to demonstrate the effective 
damping: initial energy in the trunk, dissipation in the 
branches and effective damping evaluated by the total en- 
ergy loss over one period of the trunk mode. Note that 
this model incorporates several differences with the pre- 
vious one: a very large number of modes, symmetric and 
non-symmetric modes, non-localized mass and stiffness. 



A. Model 

The model consists of three assembled beams, fig- 
Each beam has a uniform circular cross-section 



6(a) 



and is made of a linearly elastic, isotropic and homoge- 
neous material. The trunk, of length li and diameter di, 
is clamped at the base. Two symmetrical branches, each 
of length I2 and diameter c?2, are clamped at the tip of 
the trunk so that they each form an angle 0b with the 
trunk direction. As in Section |IT1 we analyze the free 




(a) 





(d) 



FIG. 6: The continuous model of a Y-shape. (a) Geometry, 
(b) Static initial condition, (c) Trunk mode, (d) Damped 
branch mode. 

vibrations of the structure. By analogy with [HI the ini- 



tial condition is an initial deformation resulting from a 
horizontal static pull on the tip of the trunk, figure |6(b)[ 
This initial energy Eq is normalized so that it is equal to 
1 when the deflection of the trunk is equal to its length 
such that X — li. To solve the equations of motion, finite- 
element computations are performed using the CASTEM 
V.3M software (Verpeaux, Charras and Millard 1988). In 
order to take into account large displacements, a step- 
by-step procedure is needed implying an update of the 
deformed configuration. Timoshenko-beam elements are 
used taking into account the rotational inertia of the sec- 
tions of the beams. Ten mesh-elements per beam were 
chosen according to a convergence test on the initial en- 
ergy repartition on modes. It appears that this number 
is sufficient to describe the full dynamics of the system. 

The first two modal shapes are given by modal anal- 
ysis and are shown in figures 6(c) and 6(d)[ We observe 



that the branch mode, figure 6(d) 



involves only branches 
bending while the trunk mode, figure [6(c)| involves both 
trunk and branch deformations. This modal localiza- 
tion and the initial condition implies that the branch 
mode has no initial energy. Note that the initial en- 
ergy is mainly in the trunk mode but can be shared 
with other modes than the branch mode. By analogy 
with Section [III two dimensionless parameters are cho- 
sen: the frequency ratio f2 = UJ2/0J1 and a mass ratio 
F = (lim2) / (hfTT-i) , where wi, 102, and toi, m2 are the 
modal angular frequencies and the modal masses of the 
trunk mode and branch mode respectively. 

As in Section |lll energy dissipation is introduced ar- 
tificially on the branch mode only by a damping matrix 
derived from the mass matrix of the finite-element model 
with a damping rate denoted ^b- The resulting effective 
damping mechanism is studied with the same definition 
of the effective damping rate, ^eff, as in Section HTbI and 
for the same reference values of the design parameters 
0b = 7r/2, ^b = 0.2, 17 = 2, and F = 0.2. 



B. Results 

The simulated dynamics of the continuous Y-shape 
yields a similar time evolution of the total energy to that 
of the lumped parameter model, figure [3] The corre- 
sponding effective damping rate is plotted in figure [7] as 
a function of the normalized initial energy Eq. As in fig- 
ureUl the effective damping rate increases with the initial 
energy level, quasi-linearly at first and reaching several 
percent for high levels of initial energy. This is a first 
indication of the robustness of the effect of branching on 
damping in a more realistic structure. As in Section |lTl 
the effects of the design parameters 0b, Cbj and F on 
the effective damping scaled by the initial energy, ^off / Eq, 
are represented in figures [51^ a)-(d) for three initial energy 
levels. Some differences appear for the continuous model 
as expected. Firstly, the branching effect on the effective 
damping is clearly maximal for larger branching angles, 
0b ~ 27r/3, rather than 7r/2, figure 8(a) Secondly, the ef- 
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Ceff 2 




FIG. 7: Effect of the initial energy level, Eo, on the effective 
damping, ^eft , of the continuous model of a Y-shape. The 
design parameters are set to ~ 7r/2, ~ 0.2, (1 = 2, and 
r = 0.2. 




FIG. 8: Effects of the design parameters on the effective 
damping scaled by the initial energy, ^es/Eo, of of the con- 
tinuous model of a Y-shape:: -Eo ~ 0.01 ( — ); Eo — 0.1 { ); 

and Eo — 1 {■ ■ ■). Unless varied, the design parameters values 
are: 0b = tt/2, £_b = 0.2, = 2, and F = 0.2. (a) Effect of the 
branching angle (j}h. (b) Effect of the branch mode damping 
^b- (c) Effect of the branch/trunk modal frequencies ratio O. 
(d) Effect of the branch/trunk modal mass ratio F. 



fective damping stabilizes or even slightly increases with 
the branch mode damping ^b, figure [8(b)| instead of de- 
creasing after ~ 0.2. Thirdly, the effective damping 
is higher at a modal frequency ratio = 3 than $7 = 2, 
figure 8(c) suggesting a richer pattern of internal reso- 



nances. Singularly, the effective damping disappears for 
n = 1, which is a striking difference with the classical 
tuned-mass damper model. Finally, the effective damp- 
ing increases with the modal mass ratio F but in a more 
complex way, figure |8(d)[ 

In conclusion, it appears that the main features of the 
damping-by-branching mechanism are still present in this 
more realistic Y-shape structure. 



IV. TWO BIOINSPIRED BRANCHED 
STRUCTURES 

Based on the results of Sections |lT] and lllli two bioin- 
spired branched structures are considered, figure O 



(a) 




(d) 





(f) 



(g) 



FIG. 9: Two bioinspired branched structures, (a) A two- 
orders ramified T-shaped structure, and its modes of interest: 
(b) Trunk mode, (c) Damped branch mode, (d) A double- 
branched Y-shaped structure, and its modes of interest: (e) 
Trunk mode, (f) Damped large branch mode, (g) Damped 
small branch mode. 



The first bioinspired structure, shown in figure 9(a)[ 
is a two-generation T-shaped structure designed so that 
4'h = 7r/2 at each branching. The ratios of branch length 
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and diameter are respectively the same between orders 
of branching, i.e. I2/I1 = h/h and ^2/0^1 = d^/d2- They 
are chosen so that the modal frequency ratio between 
the trunk mode, figure |9(b)[ and the last-order branch 
mode, figure 9(c)[ is 1:2 and so that the modal mass ra- 
tio is 0.2. A damping rate of 0.2 is introduced in this 
last-order branch mode only. The second bioinspired 
structure, shown in figure |9(d)[ consists of a double Y- 
shaped pattern with an added level of branching at 3/4 
of the height of the trunk. Both levels of branching have 
a branching angle (p^, — 27r/3 and are designed so that 
the modal frequency ratio between the trunk mode, fig- 
9(e) and the large branch mode, figure [9(f)| is 2, and 
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the modal frequency ratio between the trunk mode and 
the small branch mode, figure 9(g) is 3. A damping of 



0.2 is introduced in the two branch modes only. 

The motion is computed using the same finite-element 
method with the same initial condition, and the same 
normalization for the energy, as in Section IIIII The re- 
sulting effective damping is studied with the same damp- 
ing criterion as in Section III Bl and is plotted in figure ITUI 
as a function of the normalized initial energy. In both 




FIG. 10: Effect of the initial energy level, Eq, on the effective 
damping, ^off, of the two bioinspired structures. ( — ) Ram- 
ified T-shaped structure; ( ) Double-branched Y-shaped 

structure. 

structures, the effective damping reaches several percent 
(~ 3%) for high levels of initial energy, roughly corre- 
sponding to a third of the initial energy being dissipated 
after one period of the first mode. These results on two 
different bioinspired branched structures show that the 
damping-by-branching mechanism seems to be robust re- 
garding the branching scheme. 



V. DISCUSSION AND CONCLUSION 

At this stage, one may consider the results of the pre- 
ceding sections in relation to the proposed aim of the 



paper: to identify and characterize the elementary mech- 
anism that causes nonlinear modal energy transfer and 
damping in branched structures. In Section HIl we have 
shown that branching is the key ingredient needed to 
obtain the modal energy transfer and the resulting effec- 
tive damping that several authors had suspected. Sec- 
tions Uni and [iVj confirm that the essential features of this 
damping-by-branching are present even in more compli- 
cated branched models. Since this nonlinear mechanism 
originates in geometrical effects, the larger the amplitude 
of motion, the higher the damping. Efficient use of this 
mechanism can be made with modest requirements of 
the design parameters: a modal frequency ratio between 
the trunk mode and the branch mode near 1:2, a branch 
damping rate larger than 0.2, and the highest possible 
modal mass ratio. All these results suggest that modal 
energy transfer and the resulting damping-by-branching 
is a robust effect at large amplitudes of motion: it may 
therefore be used in a bioinspired perspective for effi- 
ciently damped slender structures. Clearly, the mecha- 
nism found here is distinct from the classical tuned-mass 
damper, hypothesized by Spatz et al. (2007) to exist in 
trees. The present damping-by-branching mechanism dif- 
fers in two essential points: (i) it is strongly amplitude- 
dependent; (ii) it is not associated with the condition of 
identical modal frequencies of the trunk mode and of the 
branch mode. 

Before generalizing our results, discussion is needed of 
some of the assumptions made to derive them. Firstly, 
all the analyses pertaining to the effective damping have 
been made on the dynamical responses to puU-and- 
release initial loading of the undamped trunk modes. 
This choice was made so that only the nonlinear geo- 
metrical effects could cause the effective damping of the 
structure. If other classical types of loading were consid- 
ered such as an initial impulse, or a harmonic or random 
forcing (Humar 2002), energy would have been given to 
all the modes of the branched structure. Although this 
would make the global energy balance more complex to 
analyze, the nonlinear geometric terms responsible for 
the modal energy transferjTjwould be still present but the 
resulting effective damping would not be simply quanti- 
fied. Secondly, we have always considered perfectly sym- 
metric and plane structures. If asymmetry between the 
branches is introduced in the model of Section|lll a linear 
coupling is introduced between the trunk and branch an- 
gles of motion so that the energy balance analysis is again 
more complex. Similarly, if three-dimensional effects are 
introduced, such as torsion or multiple 3D branching as 
in real trees, a much larger number of degrees of freedom 
is needed in order to describe the dynamics of the struc- 
ture. However, the results of Sections IIIII and IIVI show 
that complicating the modal content of the model does 
not impact the existence of the modal transfer mech- 
anism and the resulting effective damping. Moreover, 
Rodriguez et al. (2008) showed that there exist no signif- 
icant differences between the dynamics of an actual tree 
architecture and that of an idealized one. 
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More generally, the question of how branched systems 
move in a fluid environment is important in practice. In 
fact, our analysis on damping originated in the dynamics 
of trees under wind- loading (Spatz et al. 2007). In the in- 
teraction between a branched structure and flow, several 
distinct effects may be expected (Blevins 1990, Pai'doussis 
et al. 2011). Firstly, even in the absence of flow, just the 
presence of a fluid around the structure causes damp- 
ing. This damping is present in all modes, is amplitude- 
dependent, and introduces nonlinear coupling between 
modes. A flow-induced damping may also appear in 
addition. These effects may be gathered under the 
generic term of aeroelastic or hydroelastic damping. Sec- 
ondly, added mass and added stiffness effects appear 
and may alter the essential dynamical characteristics of 
the branched structure, such as frequencies and modal 
shapes. These effects are more pronounced in water. Fi- 
nally, flow may cause a large variety of loadings through 
mechanism such as turbulence excitation or wake interac- 
tions. From this list, it appears that all the key parame- 
ters involved in the mechanism of damping-by-branching 
are affected by a fluid environment: modal dampings, 
frequencies, modal shapes, and external excitations. A 
systematic analysis of these effects is clearly needed on 
the basis of the simple framework presented in this arti- 
cle. 

In the design of a slender structure, a bioinspired 
damping effect can easily be obtained according to the 
present results. The main design rule to follow is to intro- 
duce branching. Secondary rules, aiming at optimizing 



the efficiency of this damping mechanism, are to set ratio 
between the modal frequency of the branch mode and the 
trunk mode near 1:2 and the modal mass ratio as high 
as possible. Although a design rule based on the ratio of 
modal frequencies and masses is not common, it should 
be noted that the requirements are not strict, we have 
shown in Sections HIl and [IIII that damping- by-branching 
is robust and is significant for a wide range of modal fre- 
quency and mass ratios. The 1:2 rule is only indicative, 
as it is related to the original internal resonance condition 
between the branch and trunk modes. At this point, the 
concept of damping-by-branching has only been demon- 
strated to exist theoretically and numerically. Although 
it has been inspired by the observation of natural sys- 
tems, it evidently needs to be explored experimentally 
on man-made structures. 
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